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ABSTRACT 


In  1983,  continuous  acoustic  transmissions  centered  at  133  Hz  and  with 
a  resolution  of  60  ms  were  transmitted  for  five  days  from  Oahu  to  the  coast  of 
Northern  California  (4000  km  range).  A  maximum  likelihood  estimate  of  the 
change  in  acoustic  travel  time  (based  on  phase)  between  received  pulses  is  used 
to  estimate  barotropic  fluctuations.  Analysis  of  the  resulting  time  series  reveals 
resonant  oscillations  at  nontidal  frequencies  in  the  Northeast  Pacific.  Some  of  the 
periods  of  the  resonant  oscillations  are  consistent  with  theory  (Platzmzin,  Curtis, 
Hansen,  and  Slater,  1981). 

A  Wiener  filter  is  formulated  for  estimating  the  barotropic  tides  from  a  basin 
scale  tomographic  array.  Error  analyses  indicate  an  ability  to  estimate  barotropic 
currents  and  surface  displacements  with  errors  less  than  0.01  cm/s  and  1  cm, 
respectively,  over  a  large  portion  of  the  Northeast  Pacific. 
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Chapter  1 


Introduction 


1.1  Motivation 

There  are  many  reasons  to  measure  the  large  scale  barotropic  motions  of  the 
ocean.  Some  reasons  are: 

1.  Basin  scale  barotropic  motions  which  resonate  in  the  world  oceans  have 
been  theoretically  predicted  for  simplified  geometries  by  numerous  investi¬ 
gators  (eg.  Lamb,  1932,  Longuet-Higgins  and  Pond,  1970,  and  Pedlosky, 
1987).  Miller  (1986)  studied  barotropic  motions  over  variable  topography  in 
simplified  basin  shapes  by  numerically  solving  the  shallow  water  equations 
for  barotropic  vorticity  modes.  Platzman  (1979)  developed  a  finite-element 
model  of  the  world  ocean  which  allows  for  topography  and  numerically  solved 
for  the  barotropic  normal  modes  predicting  56  modes  with  periods  between 
8  and  80  hours  (Platzman,  Curtis,  Hemsen,  and  Slater,  1981).  Lacking  from 
the  studies  of  barotropic  modes  are  direct  observations.  With  ocean  acous¬ 
tic  tomography,  it  is  possible  to  measure  the  excited  modes  and  gain  insight 
into  the  forcing  mechanisms,  whether  they  be  the  wind  (Luther,  1983)  or 
the  mesoscale  eddy  field  (Miller,  1986) . 
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2.  Numerical  tidal  models  with  vaxying  degrees  of  sophistication  and  assump¬ 
tions  (eg.  Estes,  1977,  Parke  and  Hendershott,  1980,  and  Schwiderski,  1978) 
exist,  but  are  limited  to  a  tidal  height  accuracy  of  10  cm  (Schwiderski,  1978) 
and  give  no  estimates  of  the  baxotropic  currents.  Accurate  estimates  of 
surface  elevation  are  useful  in  a  number  of  applications  including  the  deter¬ 
mination  of  satellite  orbits,  mean  sea  surface  computation  (Marsh,  Martin, 
McCarthy,  and  Chovi^z,  1980),  and  satellite  measurement  of  surface  cur¬ 
rents.  With  acoustic  tomography,  tidal  height  estimates  are  improved,  and 
barotropic  tidal  current  estimates  are  a  direct  result  of  the  analysis. 

1.2  Overview 

In  Chapter  2,  we  develop  the  maximum  likelihood  estimate  of  the  tomographic 
travel  time  fluctuations  due  to  barotropic  motions.  We  apply  this  technique  to 
five  days  of  data  from  the  jw:oustic  tomography  experiment  of  1983  and  perform 
spectral  analysis  on  the  resulting  travel  time  series. 

In  Chapter  3,  we  formulate  a  Wiener  filter  to  estimate  the  barotropic  tides. 
Error  maps  based  on  synthetic  data  are  shown. 
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Chapter  2 


Measurement  of  Travel  Time  Changes  Due  to 
Barotropic  Motions 


The  acoustic  tomography  experiment  in  1983  (abbreviated  to  AT83)  was  con¬ 
ducted  by  Spiesberger,  Metzger,  and  Birdsall  and  is  the  subject  of  two  papers  in 
preparation  (Spiesberger,  Bushong,  Metzger,  ajid  Birdsall,  1987,  and  Spiesberger, 
Metzger,  and  Birdsall,  1987).  V-^e  present  a  short  description  of  the  AT83  exper¬ 
iment,  a  method  for  measuring  the  tomographic  travel  time  changes  caused  by 
barotropic  motions,  and  the  results  of  this  amalysis  on  the  AT83  data. 


2.1  Geometry  and  Signal  Processing  of  the  AT83 
Experiment 

The  AT83  experiment  consisted  of  approximately  one  week  of  acoustic  trans¬ 
missions  between  a  bottom- mounted  source  at  183  meter  depth  on  the  coast  of 
Oahu  and  a  bottom-mounted  receiver  at  mid-depth  off  the  coast  of  California 
(Figure  2.1).  The  distance  W2is  about  4000  km.  The  transmitted  signal  had  a 
center  frequency  of  133  Hz  and  a  bandwidth  of  about  17  Hz.  The  source  was  not 
able  to  transmit  a  pulse  which  could  be  received  above  the  noise  at  the  receiver. 

I 
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Figure  2.1:  The  acoustic  tomography  experiment  in  1983  consisted  of 
a  source  (S)  located  neair  Kaneohe  Bay,  Oahu,  and  a  bottom-mounted 
receiver  (R)  whose  approximate  position  is  shown.  The  distance  be¬ 
tween  the  instruments  was  about  4000  km. 


Instead,  the  133  Hz  signal  was  phase-modulated  every  8  cycles  using  a  511  digit 
TnaviTna.!  shift  register  sequence  having  a  period  of  511(8/133)  =  30.7  seconds. 
Continuous  transmissions  were  made  between  yeardays  333  (29  November)  zmd 
339  (5  December)  in  1983  except  for  two  gaps,  one  of  2  hours  duration  and  the 
other  of  5  hours  duration. 

The  received  signal  was  complex  demodulated  and  correlated  with  a  phase 
modulation  replica  (phase-only  filtering).  Replica  correlation  produced  an  output 
equivalent  to  the  signal  which  would  be  received  if  a  60  ms  pulse  was  transmitted 
every  30.7  seconds.  This  correlation  also  gave  a  signal-to-noise  ratio  increase  of 
lOloKio^ll  =  27  dB  along  each  ray  path.  Four  successive  pulse  receptions  were 
averaged  together  to  increase  the  signal-to-noise  ratio  by  an  additionad  10  log^Q  4  = 
6  dB  and  resulted  in  multipath  records  at  two  minute  intervals.  The  final  record 
was  sampled  at  about  15  ms  intervals  (two  cycles  of  carrier).  A  running  average 
across  fom  samples  increased  the  signad-to-noise  ratio  by  10  log^g  4  =  6  dB  with 
the  loss  of  some  resolution.  These  processing  steps  are  identical  to  those  used 
before  (Spindel,  1979,  and  Spiesberger,  Spindel,  amd  Metzger,  1980)  except  for 
the  specific  frequencies,  averaging  intervals,  amd  shift  registers  used. 

Multipaths  arrived  over  a  four  second  interval.  Adjacent  acoustic  records  do 
not  look  similar  because  of  the  fluctuations  of  internal  waves  (Spiesberger,  Metzger, 
and  Birdsadl,  1987).  Incoherent  averages  over  a  period  equad  to  the  lowest  internad 
wave  period  (about  one  cycle  per  day  at  this  latitude)  filter  out  much  of  the  short 
term  acotistic  variability  (Spiesberger,  Spindel,  and  Metzger,  1980).  Figure  2.2 
shows  the  result  of  incoherently  averaging  703  acoustic  records  (24  hours  of  data) 
on  each  of  four  days.  There  are  five  stable  ray  airrivads  (A-E). 

The  time  base  of  the  signal  generator  on  Oahu  was  maintained  by  a  Rubidium 
oscillator  and  a  GOES  satellite  receiver.  The  time  at  which  the  signal  left  the  signal 
generator  was  measured  to  an  accuracy  of  1  ms.  The  signal  travelled  through 
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Figure  2.2:  Incoherent  averages  of  received  pulses  over  24  hours  (703 
pulsed  transmitted  at  intervals  of  about  2  minutes).  The  60  ms  pulses 
transmitted  at  the  Kaneohe  source  are  spread  by  multipaths  and  ar¬ 
rive  at  the  Northern  California  coast  over  an  interval  of  about  four 
seconds.  Five  stable  ray  arrivals  (A-E)  are  indicated. 


13 


a  power  amplifier,  a  5  mile  length  of  cable,  zmd  through  the  acoustic  source. 
The  delay  through  the  acoustic  source  is  known  to  zm  accuracy  of  about  20  ms. 
Although  the  transmission  time  is  only  known  to  20  ms,  the  change  in  transmission 
time  was  known  to  about  10  /xs  during  the  experiment. 

The  time  base  at  the  receiver  was  also  maintained  by  a  Rubidium  oscillator 
and  a  GOES  satellite  receiver.  The  accuracy  of  the  receiver  clock  was  about  2  ms. 
The  drift  of  this  clock  was  less  than  10  fxs  during  the  experiment. 

Since  the  source  and  the  receiver  were  bottom-mounted,  acoustic  travel  time 
changes  exceeding  about  14  /xs  were  due  to  ocean  fluctuations  or  noise. 


2.2  The  Effect  of  Barotropic  Motions  on 
Tomographic  Travel  Times 


The  change  in  acoxistic  travel  time  along  a  ray  path  To  due  to  fluctuations  in 
temperature,  6$,  sea  surface  displacement,  rjo,  and  currents,  u,  is  given  by. 


«r(t)  =  -a/  -  kf  (2.1) 

Jro  cs(x,t)  ho  cl{x,t)  Jto  c5(i,f) 

(Munk  and  Wunsch,  1978,  and  Munk,  et  al.,  1981).  The  position  along  the  ray 
path  is  X,  and  the  geophysical  time  is  t.  An  increment  along  the  ray  path  is  ds  and 
has  a  direction  given  by  s.  The  reference  sound  speed  is  denoted  by  cq.  Typical 
values  for  k  and  a  are  0.017  s"^  and  4.785 


Barotropic  fluctuations  whose  horizontal  scale  is  much  greater  than  a  ray  cy¬ 
cle  distance  (about  50  km)  will  siflect  all  ray  travel  times  about  the  same.  The 
barotropic  tides  have  large  scales  (order  1000  km).  Typical  values  are. 


f/o  ~  1  meter 
u  ~  0.02  m/s. 
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The  associated  isotherm  displacement  causes  a  temperature  of  about 


Sd  -  0.005“ C 

in  the  main  thermocline.  The  ray  will  spend  about  one  quarter  of  its  time  in  the 
displaced  thermocline.  For  a  range  of  4000  km,  the  magnitudes  of  the  three  terms 
in  equation  (2.1)  are  10  ms,  30  ms,  and  36  ms,  respectively.  In  the  next  section 
we  demonstrate  that  these  perturbations  are  significzmt  In  comparison  with  the 
imcertainty  of  our  travel  time  measurements. 


2.3  Measurement  of  Travel  Time  Change  Using 
Phase 

We  seek  a  method  of  estimating  the  average  travel  time  change  between  records 
in  order  to  search  for  large  scale  barotropic  processes.  Travel  time  measurements 
based  on  phase  are  usually  more  accurate  than  those  based  on  amplitude. 

Acoustic  travel  time  changes  appear  as  changes  in  the  phase  of  the  demodulated 
complex  signal;  the  travel  time  change,  r,  being  linearly  related  to  the  phase,  (f>, 
through  0  =  cjr  where  w  is  the  carrier  frequency.  Measured  phase  is  a  modulo  1 
cycle  process.  Ambiguities  may  result  when  extending  the  phase  past  one  cycle 
(Dyson,  Munk,  and  Zetler,  1976). 

The  l/e  decorrelation  time  of  the  complex  cross  correlation  for  our  records  is 

4.3  minutes  (Figure  2.3).  We  can  use  the  high  correlation  of  adjacent  2u:oustic 
signals  (separated  by  two  minutes)  to  estimate  the  average  chzmge  in  phase.  The 
two  adjacent  complex  records  can  be  written  as. 


=  fi{t)  exp  (j(w«t  +  (f>i  +  fl(t))] 

(2.2) 

r2(t) 

=  /2(0exp[j(wet  +  <^2  +  fl(t))], 

(2.3) 

15 

Correlation 


Correlation  Lag  (Min) 


Figure  2.3:  Peak  cross  correlation  magnitudes  versus  transmission 
time  separation  indicate  a  1/e  decorrelation  time  of  4.3  minutes.  Each 
point  is  an  average  zurross  ten  realizations. 
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where  fi{t)  and  f2{t)  are  non-negative  valued  functions,  (f>2  ~  4>i  is  the  average 
change  in  phase,  and  di(t)  and  d2{t)  account  for  the  other  phase  fluctuations.  We 
wish  to  estimate  the  change  in  the  average  phase,  A«^2  =  every  two  minutes. 

The  sufficient  statistics  X  and  K,  defined  as. 


where  *  indicates  the  complex  conjugate,  lead  to  the  maximum  likelihood  (matched 
filter)  estimate  of  the  phase  change, 

A(f>  =  arctan  (2.6) 

(Van  Trees,  1968).  The  upper  limit  of  integration  of  equations  (2.4)  and  (2.5),  St, 
is  chosen  to  include  only  the  incoming  signals.  The  unwrapped  or  extended  phase 
is  then, 

<t>n  =  (2.7) 

where  the  subscript  n  indicates  the  record  number,  and  an  initial  phase,  of 
zero  is  assumed.  We  call  this  matched  filter  processing  because  the  filtering  of  the 
second  sigpial  is  matched  to  the  previous  signal. 

Previous  continuous  wave  (CW)  experiments  have  estimated  travel  time  fluctu¬ 
ations  from  the  means  of  the  real  2md  imaginary  parts  of  the  complex  demodulated 
signal.  If  x(t)  and  y(t)  are  the  real  and  imaginary  parts  of  the  complex  demod¬ 
ulated  signal,  respectively,  aind  fit  is  the  averaging  period,  the  estimates  of  the 
means  are. 


fii 


with  a  phase  estimate  of 


4>{t)  =  axctan 


(2.10) 


Difficulty  arises  in  extending  the  modulo  1  cycle  phase  when  the  phase  chamges 
are  large. 

The  first  two  days  of  the  AT83  experiment  data  were  processed  using  both 
methods  (Figure  2.4).  CW  processing  was  performed  by  calculating  X(f)  zuid 
y’(t)  according  to  equations  (2.8)  and  (2.9).  The  limit  of  integration,  6t,  was  four 
seconds  to  include  only  the  signal.  Phase  was  extended  according  to, 


\<f>n  -  1/2(^„_1  +  ^n-:)|  <  1/2 cycle. 
Phase  kinks  were  then  removed  according  to. 


(2.11) 


-  l/2(<^n+i  +  (^n-i)l  <  1/2  cycle  (2.12) 

(Dyson,  Munk,  and  Zetler,  1976).  The  long-period  divergence  of  the  estimators  is 
caused  by  the  ambiguity  of  extending  the  CW  record  after  large  phase  changes. 
Figure  2.5  shows  histograms  of  CW  and  matched  filter  phase  chemge  over  two 
minutes  for  five  days  of  data.  The  large  phase  chzmges  in  the  CW  record  can 
be  attributed  to  the  integration  processes  of  equations  (2.8)  and  (2.9)  which  are 
sensitive  to  any  multipath  fading. 


2.4  Observations  of  Resonant  Oscillations  in  the 
Northeast  Pacific 

We  use  spectral  analysis  to  investigate  the  estimates  of  travel  time  (Figures  2.6 
and  2.7)  over  the  first  five  days  of  the  AT83  experiment.  All  spectral  estimates 
are  one-sided.  A  line  fit  was  used  to  detrend  the  travel  time  series.  The  trend 
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CYCLES 


1983 


Figure  2.4:  Two  estimators  of  the  acoustic  travel  time  change  (one 
phase  cycle  is  about  7  ms)  based  on  phase  from  two  day  of  data.  The 
jagged  curve  is  the  estimate  obtained  if  the  source  transmitted  a  CW 
signal  at  133  Hz  where  phase  is  unwrapped  according  to  equations 
(2.11)  and  (2.12).  The  smooth  curve  is  the  estimate  obtained  with 
matched  filter  processing.  The  long-period  divergence  of  the  estima¬ 
tors  is  caused  by  the  ambiguity  of  unwrapping  the  CW  record  after 
large  phase  changes. 
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Figure  2.5:  Histograms  of  the  measured  phase  change  in  two  minutes 
over  five  days  of  data  (3486  records)  using: 

Top:  Matched  filter  processing 
Bottom:  CW  processing. 
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Figure  2.6:  The  matched  filter  processing  used  to  estimate  the  acous¬ 
tic  travel  time  change  over  4000  km  from  the  Kaneohe  source  near 
Oahu  to  the  Northern  California  coast. 
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Figure  2.7:  Acoustic  travel  time  change  of  Figure  2.6  detrended  by  a 
line  fit.  Over  the  five  day  series,  the  trend  accounts  for  a  travel  time 
increase  of  38  ms  corresponding  to  a  cooling  of  the  Northeast  Pacific. 
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1983 


Spectral  Technique 

Peak  Number 

Period  (hr) 

Varizmee  (cycles*) 

Periodogram 

1 

1.8 

Autoregressive 

1 

21.2 

1.7 

Prony 

1 

21.4 

1.3 

Periodogram 

2 

15.8 

1.6 

Autoregressive 

2 

15.5 

2.4 

Prony 

2 

15.7 

3.7 

Periodogram 

3 

12.2 

5.0 

Autoregressive 

3 

12.3 

6.0 

Prony 

3 

12.3 

4.3 

Periodogram 

4 

9.8 

— 

Autoregressive 

4 

9.6 

— 

Prony 

4 

9.4 

1.5 

Table  2.1:  A  summary  of  the  prominent  peaks  from  three  spectral  estimates  of 
the  travel  time  series  of  Figure  2.7 

accounted  for  a  38  ms  increase  in  acoustic  travel  time  corresponding  to  a  cooling 
of  the  Northeast  Pacific. 

We  use  three  spectral  estimators  tc  measure  the  frequencies  of  variation  of  the 
travel  time  series;  a  periodogram  (Figure  2.8),  the  autoregressive  method  (Figure 
2.9),  and  Prony’s  method  (Figure  2.10).  The  periodogram  is  a  classical  spectral 
technique  which  is  robust  for  long,  wide-sense  stationary  series,  but  is  limited  by 
sidelobes  zmd  resolution  for  short  series.  Autoregressive  and  Prony’s  method  are 
more  suitable  for  short  series  because  they  have  high  resolution  without  sidelobes 
(Kay  and  Marple,1981,  and  Marple,  1987). 

We  discuss  four  peaks  which  all  three  spectral  estimators  exhibit  (Table  2.1). 
Four  observations  are  made: 

1.  The  diurnal  tide  peak  which  is  expected  at  a  period  of  about  24  hours  is  not 
present.  A  peak  at  a  period  near  21.2  hours  exists. 


2.  A  peak  occurs  at  a  period  of  15.7  hours. 


Fowcn*  Spectral  Density 


Figure  2.8:  Periodograru  estimate  of  the  power  spectral  density  (cycles* /cph) 
of  the  detrended  phase  series  (Figure  2.7).  The  three  major  peaks  oc¬ 
cur  at  periods  of  20.7,  15.8,  and  12.2  hours. 
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Figure  2.9:  Autoregressive  method  estimates  of  the  power  spectral 
density  (cycles* /cph)  of  the  detrended  phase  series  (Figure  2.7). 

Top:  Sample  lags  of  30.7  min  and  a  model  order  of  75. 

Bottom:  Sample  lags  of  6.1  min  and  a  model  order  of  400. 

In  both  estimates,  the  major  peaks  occur  at  periods  of  21.2,  15.5,  amd 
12.3  hours. 
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Figure  2.10:  Prony’s  method  estimates  of  the  energy  spectral  density 
(cycles* /cph*)  of  the  detrended  phase  series  (Figure  2.7). 

Top:  Sample  lags  of  30.7  min  and  a  model  order  of  75. 

Bottom:  Sample  lags  of  6.1  min  and  a  model  order  of  400. 

In  both  estimates,  the  major  peaks  occur  at  periods  of  21.4,  15.7,  and 
12.3  hours. 


26 


Mode  Niimber  Period  (hr) _ Mode  Description 


16 

28.7 

Pacific  1/2  wave;  Antarctic  Kelvin  wave  1 

19 

21.2 

Pacific  2/2  wave 

25 

15.5 

Pacific  North  American  Kelvin  wave 

26 

15.2 

Antarctic  Kelvin  wave  2 

32 

12.5 

Arabian  Sea;  North  Pacific 

35 

11.6 

Equatorial  Pacific  transverse  3/2  wave 

46 

9.2 

Equatorial  Pacific  transverse  4/2  wave 

Table  2.2:  Predicted  barotropic  normal  modes  of  the  Northeast  Pacific  with  peri¬ 
ods  between  9  and  30  hours  (from  Platzman,  Curtis,  Hzmsen,  zmd  Slater,  1981). 

3.  The  semi-diurnal  tide  peak  which  is  expected  near  a  period  of  12  hours  is 
present. 

4.  A  smaller,  but  significant,  peak  appears  at  a  period  of  9.2  hours. 

In  seeking  to  explain  these  results,  we  cite  the  work  of  Platzman  (Platzman, 
1979, 2md  Platzman,  Curtis,  Hansen,  and  Slater,  1981).  He  developed  a  barotropic 
finite-element  model  of  the  world  ocean  which  allowed  for  topography  and  numer- 
icadly  solved  for  the  barotropic  normal  modes,  finding  56  modes  with  periods  be¬ 
tween  8  and  80  hours.  Of  these  56  modes,  41  have  periods  between  8  and  30  hours 
and  are  primarily  gravity  modes,  that  is,  gravity  provides  the  restoring  force.  A 
listing  of  some  of  Platzman ’s  normal  modes  with  periods  between  9  and  25  hours 
is  given  in  Table  2.2. 

The  spectral  peaks  at  periods  of  21.2, 15.5,  and  9.2  hours  correspond  to  Platz- 
man’s  modes.  The  phase  series  length  is  insufficient  to  resolve  both  the  15.2  and 
15.5  hour  period  components.  Other  significant  peaks  appear  in  all  three  spectral 
estimators  (periods  of  7.5, 6.4,  and  5.3  hours)  but  have  periods  less  than  those  com¬ 
puted  by  Platzman. 
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Chapter  3 


Estimation  of  Barotropic  Tides 

We  formulate  an  inversion  filter  to  estimate  the  barotropic  tides.  The  periodic 
barotropic  motions  are  required  to  satisfy  conservation  of  mass  and  are  specified 
by  fo\ir  current  models.  The  earth  is  assumed  to  be  spherical,  and  vertical  dis¬ 
placements  are  relative  to  the  earth’s  surface.  Root  mean  square  estimation  errors 
are  explored  for  a  source-receiver  array  similar  to  that  of  the  aicoustic  tomography 
experiment  of  1986-87. 

3.1  Forward  Problem 

The  change  in  acoustic  travel  time  along  a  ray  path  To  due  to  isotherm  dis¬ 
placement,  rj,  surface  displacement,  rjo,  and  currents,  u,  follows  from  equation 
(2.1)  and  is  given  by, 

6T{t)  =  •  3  -I-  ds,  (3.1) 

where  dcp/dz  is  the  potential  sound  speed  gradient  (Mimk,  et  al.,1981)  and  z  = 
ro  —  r  where  ro  is  the  earth’s  radius.  A  representative  value  of  soimd  speed,  c,  has 
been  substituted  for  the  position  dependent  sound  speed. 


r 


The  barotropic  current  as  a  function  of  colatitude,  6,  longitude  (east),  (f>,  fre¬ 
quency,  (j,  and  time,  t,  is  taken  to  be, 

u(^,^,w,t)  =  (ai(d,0,u;)costi;t  +  Oj(d,0,a;)  sinwt]  d 

+  [6i(d,^,w)  coswt  +  b2{0f(t>,u)  sinwt]  $.  (3.2) 

Surface  elevation  and  currents  are  related  through  the  continuity  of  mass  equation, 

^2^  =  -V*  .  (k(i)ff(x,()) ,  (3.3) 

where  h  is  the  ocean  depth  (positive)  and  Vj^  is  the  horizontal  component  of 
the  gradient.  Substituting  equation  (3.2)  into  equation  (3.3)  and  integrating  with 
respect  to  time  gives. 


— ^  \hai  cot  0  +  half  +  hbxA  esc  B  +  htai  +  h^bi  esc  B]  sinwt 
ru 

[Aaj  cot  B  +  ha2f  +  hb^^  esc  B  +  hfOt  +  h^bz  esc  B\  cos  wt,  (3.4) 

where  the  subscripts  B  and  <l>  indicate  the  derivative  with  respect  to  colatitude  and 
longitude. 

To  relate  the  depth  dependent  vertical  displacement  to  currents,  we  start  with 
the  continuity  equation  in  terms  of  spherical  component  velocities,  u,  v,  and  w. 


dw  2w  1 
dr  r  rsintf 


(w  sin  B]f  + 


rsin^ 


=  0. 


From  equation  (3.2), 


(3.5) 


V  =  ai  cos  u/t  ■+■  at  sin  uft  (3.6) 

u  =  6iCosu;t -H  63sinb;t.  (3.7) 

Substituting  these  expressions  Into  equation  (3.5)  and  solving  the  differential  equa¬ 
tion  with  the  botmdary  condition  of  zero  normal  flow  at  the  bottom  gives  the 
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verticzd  velocity  as  a  function  of  position  and  time.  Integrating  with  respect  to 
time  gives  the  depth  dependent  vertical  displacement  as  a  fimction  of  position, 
frequency,  and  time, 

r}{d,<f>,z,u,t)  = 

—  z){ai  cot  d  +  tti#  +  bi^  CSC  6)  +  hgai  +  h^bi  esc  0\  sin  wis 
+—  [{h  —  z){a2Cotd  +  a^t  +  62^  esc  d)  +  htoi^  +  ^^62  esc  cos  wt  (3.8) 

TUI 

where  z  is  positive  downwards  with  z  =  0  at  the  mean  surface.  The  depth  indepen¬ 
dent  portion  of  equation  (3.8)  is  the  surface  displacement,  rfo{x,t).  Substituting 
equation  (3.8)  into  equation  (3.1)  gives  the  frequency  and  time  dependent  travel 
time  fluctuation, 

6T{t,uj)  =  {k  -  r?o(x,t)ds  -  i /^^u(g,0  -  Ids 

f  ( )  (fli  cot  9  +  ait  +  b^,  esc  0)  sinu/t  ds 
c*ruj  JTo  \  az  j 


c^rw  JTo  \ 


acp(x. 


(02  cot  0  +  02»  +  62^  CSC  $)  COS  ut  ds.  (3.9) 


Evaluation  of  equation  (3.9)  requires  path  integration  along  single  rays.  We 
remove  this  requirement  by  approximating  equation  (3.9)  with  sin  integr2d  in  the 
horizontal  only  by  assuming  that  the  path  integrals  and  horizontal  integrzils  czm 
be  related  through  constants  of  proportionality  so  that, 

ST{t,uj)  =  J  ijo{x,t)dl- ^  j  u{x,t)-sdl 

1  ^  A  f 

-1-:= —  z-:^]  I  (ai  cot  9  +  ail  +  bi4  CSC  B)  sin utdl 
c^ruj  \  az  I  J 


(‘^’^1  y  {®jC0t5 -h  02*  +  62^csc5)  coswtd/,  (3.10) 


where. 


=  /ro<^ 

^  -  Jdl 


(3.11) 
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(3.12) 


dz  S  dl 

dz  ~  Jdl 


(3.13) 


and  dl  is  the  increment  in  range  along  the  transmission  path.  For  notational 

convenience,  we  define,  _ 

Cl  =  fcTf  -  ^  (3.14) 

cj  =  (3.15) 

and  6{x)  equals  the  bearing  to  the  receiver  at  any  position,  x,  adong  the  transmis¬ 
sion  path.  We  substitute  equations  (3.2)  and  (3.4)  into  equation  (3.10)  to  get  the 
travel  time  fiuctuation  as  a  function  of  time  and  frequency, 


6Ti{t,i4j)  =  e,(u;)  coswt  +  /,(u>)  sinwt,  (3.16) 


where, 

e,-(ci;)  =  — y{(cih  +  C2)(oj  cot 0  +  Uji  +  CSC 0) -|- Ci(/i|02  +  CSC fl)}dl,- 
J  (oicostf  -  6isin^)dlj  (3.17) 


/i(w)  = 


+  C2)(aiCotd  +  ai»  +  614 esc 0)  +  Ci(hfai  +  h^bxczc6)}dli 

c*ru  J 

j  {at  coa  6  —  bi  am  6)dli,  (3.18) 


and  the  subscript  %  indicates  the  source-receiver  pair.  In  this  formulation  of  the 
forward  problem,  the  data  are  e,-(u;)  and  and  are  linearly  related  to  the 
models,  ai(d,0,(j),  02(0,0,0;),  6i(0,0,u;),  and  b2{9,<^,u)  as  in  equations  (3.17) 
and  (3.18). 
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3.2  Inverse  Problem 

We  have  chosen  to  use  a  Wiener  filter  to  estimate  the  models  (Cornuelle,  1983, 
The  Ocean  Tomography  Group,  1982).  With  N  transmission  paths,  the  2N  x  1 
data  vector, 

d(w)  =  [ei(w)  /i(u;)  e2(w)  /2(w)  •  •  •  /n(w)|^  ,  (3.19) 

is  formed  from  the  travel  time  coefficients  of  equation  (3.16). 

The  2N  X  2N  data-data  covariance  matrix, 

D'{u)=E  [d(w)c?^(w)]  ,  (3.20) 

and  the  2N  x  1  model-data  covariance  vectors, 

)3ai(5,^,w)  =  E{dai{e,<j>,ui)) 

/?!ia(^,^iw)  =  E(da2(&,<^,u)) 

=  E{db^{9,<i>,u)) 

=  E{db2{9,(l>,u,)),  (3.21) 

can  be  computed  from  equations  (3.17)  and  (3.18)  given  the  auto-correlation  and 
cross-correlation  functions  of  ai,  az,  6i,  zmd  6}.  Each  element  of  the  data-data  co- 
variance  matrix  requires  a  double  integral  over  two  transmission  paths  while  each 
element  of  the  model-data  covarizmce  vectors  requires  an  integrsJ  over  one  trans¬ 
mission  path.  The  computed  datardata  covau'iance  matrix,  IP'(u;)  is  augmented  by 
a  noise  covariance  matrix,  to  account  for  noise  in  the  data  vector,  dy  and  to 

account  for  approximations  the  forward  problem.  The  total  datardata  covariance 
matrix,  D{u;),  is  given  by, 

=  i?*(w)  +  Z?^(u;).  (3.22) 
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The  discrete  Wiener  filter  or  stochastic  inverse  estimate  of  a\{0,<i>,ui)  is  given 

by, 

(3.23) 

and  has  a  mean  square  estimation  error  of 
B  [(<ii(», *,!./)  -  <fi(#,^,u))’]  =  E  [(a,(»,^,w))’]  - 

(3.24) 

(Kailath,  1981,  and  Aki  and  Richards,  1980).  Substituting  the  appropriate  model- 
data  covariance  vector  gives  similar  expressions  for  03,61,  and  63. 

Estimates  of  the  surface  displacement,  rjo,  follow  directly  from  equation  (3.4) 
and  are  given  by, 


rfQ(tf,^,w,t)  = 

1  r  ^  ^ 

-  —  [hdi  cot  0  +  halt  +  kbit  ^  sinwt 

cot  0  +  ha\t  +  hb2t  ^  (3*25) 

The  mean  square  error  of  this  estimate  is 

E  [(t7o  -  »fo)’]  =  E  [q^]  -  E  [rioifo] ,  (3.26) 

where  i?[^o]  is  the  a-priori  variance  of  the  surface  displacement.  The  reduction  in 
variance,  .^[norfo],  follows  from  equations  (3.4),  (3.21),  and  (3.25)  and  is  given  by, 

^l^o'fo)  =  ^^^5^(w)2?”^(w)a(w),  (3.27) 

where 

o(w)  =  {ht  +  /icottf)/3a,(w)  +  +  htcac00i^{(j)  +  hcsc0-~^^~.  (3.28) 

09  0(p 
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3.3  Error  Analysis 


The  aicoustic  tomography  experiment  of  1986-87  (AT86-87)  consisted  of  three 
deep  ocean  tethered  sources,  one  source  near  Kaneohe  Bay,  Oahu,  2ind  about  ten 
bottom  mounted  receivers.  The  area  covered  by  the  zirray  was  about  4000  km  x 
3000  km.  To  explore  the  ability  of  the  inversion  filter  to  account  for  the  a-priori 
variances  of  the  currents  and  surface  displacements,  we  choose  a  source-receiver 
array  with  four  sources  that  correspond  to  those  of  the  AT86-87  experiment  but 
with  only  five  receivers  (Figure  3.1).  The  auea  covered  by  this  array  is  similar  to 
that  of  the  AT86-87  experiment,  but  with  less  resolution. 

We  assume  that  the  current  components  are  statistically  independent  and  that 
the  autocorrelation  function  of  e2u:h  current  component  is  Gaussian  ajid  a  function 
only  of  the  distance,  X,  between  the  two  geographic  points.  That  is, 

a*  exp 

a*  exp 

2  (  ^ 
a*exp(^-  — 

<7*exp^-^^,  (3.29) 

and. 


•S' [tti(^i»^i»w)ai(d2,(^a,u;)]  = 

^  (aj(^i,  0i,w)a2(d2>02>w)]  = 


E[aia2]  —  ^[<ii6ij  —  £^[0x62]  —  —  -^[<12^2)  ~  ^[6162  ~  0>  (3.30) 


where  a*  is  the  variance  of  ax,  a2,  61,  amd  62,  and  L  is  the  1/e  decorrelation  length. 

A  decorrelation  length  of  1500  km  is  assumed,  and  adl  error  analyses  are  com¬ 
puted  for  a  frequency  equal  to  that  of  the  M2  tidal  component.  The  a-priori 
variance  on  eau:h  current  model  is  tadcen  to  1.0  cm^/s^  to  give  an  a-priori  error 
of  the  current  magnitude  of  2.0  cm/s,  typical  of  the  magnitude  of  the  M2  tidal 
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Figure  3.1:  The  acoustic  tomography  experiment  of  198fr-87  consisted 
of  three  deep  ocean  tethered  sources  (S1-S3),  one  source  (S4)  near 
Kaneohe  Bay,  Oahu,  and  about  ten  receivers,  five  of  which  (R)  are 
shown  at  approximate  positions  along  the  West  coast  of  the  United 
States.  The  area  covered  by  this  array  was  about  4000  km  x  3000  km. 
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current  (Munk,  Snodgrass,  and  Wimbush,  1970).  Depths  and  bottom  slopes  are 
computed  from  Schwiderski’s  “hydrodynamical  ocean  bathymetry”  (Schwiderski, 
1978). 

The  Wiener  filter  is  able  to  estimate  the  barotropic  current  models,  Ui,  02, 
61,  and  621  accurately  over  large  areas  (Figures  3.2  and  3.3).  Currents  in  areas 
approximately  1000  km  x  1000  km  can  be  mapped  with  an  error  less  them  0.01 
cm/s. 

The  a-priori  variance  of  the  surface  displacement  can  be  evaluated  using  equa¬ 
tion  (3.4)  and  the  assumptions  of  equations  (3.30)  and  (3.31)  to  give, 

j  r  >4  L  2  * 

^bo]  =  -^  {ht  +  hcotdy  +  hlcsc^0+^  .  (3.31) 

As  evident  from  equation  (3.31),  the  natural  {a-priori)  variability  of  the  surface 
displacement  (Figure  3.4)  is  closely  related  to  the  bottom  topography  (Figure 

3.5) .  The  Wiener  filter  is  able  to  estimate  suiface  displzw:ements  over  large  areas 
(approximately  2500  km  x  2000  km)  to  an  accuracy  greater  than  1  cm  (Figure 

3.6) . 

The  inversion  filter  is  robust  in  the  presence  of  noise.  A  travel  time  uncertainty 
of  10  ms  has  minimal  effect  on  the  current  and  surface  displacement  estimation 
errors  (Figures  3.7,  3.8,  and  3.9). 
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Figure  3.2:  Root  mean  square  estimation  error  (cm/s)  of  the  north/south 
barotropic  coirrent  models  (di  and  dj)  for  a  Gaussian  decorrelation 
length  of  1500  km  and  no  data  noise.  An  a  —  priori  variance  of  1.0 
cm^/s^  is  assumed. 
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Figure  3.3:  Root  mean  square  estimation  error  (cm/s)  of  the -at/west 
bzirotropic  current  models  (bi  and  62)  for  a  Gaussian  decorrelation 
length  of  1500  km  and  no  data  noise.  An  a  —  priori  variance  of  1.0 
cm^/s’  is  assumed. 
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Figiire  3.4:  Natural  variability  (m)  of  the  surface  displacement  (rjo)  in 
the  absence  of  any  measurements  for  a  Gaussian  decorrelation  length 
of  1500  km.  An  a  -  priori  variance  of  l.O  cm*/s*  for  each  current 
model  (ai,a2,6i,and  b})  is  assumed. 
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Figiire  3.5:  Ocean  depths  (m)  for  the  area  covered  by  the  source- 
receiver  array  of  Figure  3.1.  Depth  contours  are  at  intervals  of  200 
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Figure  3.6:  Root  mean  square  estimation  error  (cm/s)  of  the  surface 
displacement  (r/o)  for  a  Gaussian  decorrelation  length  of  1500  km  and 
no  data  noise.  An  o  -  pnon  variance  of  1.0  cm*/s*  for  each  current 
model  (ai,as,6i,and  bj)  is  assumed. 
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Figxxre  3.7:  Root  mean  square  estimation  error  (rm/s)  of  the  north/south 
barotropic  current  models  (di  and  d^)  for  a  Gaussian  decorrelation 
length  of  1500  km  and  a  data  noise  level  of  10  ms.  An  a  —  priori 
variance  of  1.0  cm^/s^  is  assumed. 
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Figure  3.8:  Root  mean  square  estimation  error  (cm/s)  of  the  east /west 
barotropic  current  models  (6i  and  for  a  Gaussian  decorrelation 
length  of  1500  km  and  a  data  noise  level  of  10  ms.  An  a  —  priori 
variance  of  1.0  cm^/s^  is  assumed. 
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Figure  3.9:  Root  mean  square  estimation  error  (cm/s)  of  the  surface 
displacement  {rfo)  for  a  Gaussian  decorrelation  length  of  1500  km  and 
a  data  noise  level  of  10  ms.  An  a  —  priori  variance  of  1.0  cm*/s*  for 
each  current  model  (ai, 02,61,  and  b2)  is  assumed. 
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3.4  Conclusions 


Movement  of  the  tethered  sources  results  in  tomographic  travel  time  changes 
which  must  be  taken  into  account  when  estimating  the  data,  £{((<;)  and  fi{u)) 
(Comuelle,  1985). 

A  lack  of  trzmsmissions  for  a  period  exceeding  the  decorrelation  time  of  the 
received  pulses  causes  an  ambiguity  in  ^tending  the  travel  time  fluctuation  series. 
The  rate  of  change  of  the  travel  time  is  not  ambiguous  across  transmission  gaps. 
Since  the  barotropic  tides  are  periodic,  the  data,  stnd  can  be  computed 
from  the  differentiated  travel  time  series. 

The  time  for  sound  to  travel  4000  km  is  about  44  minutes.  Neglect  of  the  change 
in  currents  and  surface  displacement  over  this  interval  introduces  significant  errors 
in  the  estimates  (the  ctirrents  and  stirface  displacement  of  a  semi-diurnal  tide  can 
change  by  as  much  as  37  percent  of  their  amplitude).  Correction  is  made  by 
including  a  time  delay  term  (a  function  of  the  transmission  path  and  geographic 
position)  in  the  forward  problem. 

The  inversion  filter  presented  is  relatively  unconstrained,  requiring  only  con¬ 
servation  of  mass  and  assuming  statistical  independence  of  the  current  models. 
Application  of  the  momentum  equation  halves  the  number  of  independent  models 
to  estimate  and  reduces  the  estimation  errors  with  a  cost  of  additional  compu¬ 
tation.  The  estimation  errors  will  also  be  improved  by  increasing  the  number  of 
receivers  (there  were  approximately  ten  in  the  AT86-87  experiment)  (Bushong, 
Spiesberger,  1987).  Bottom  pressure  gauge  data  can  also  be  included  giving  a 
reduction  in  error. 

The  computations  required  were  done  in  a  personal  computer  on  a  Definicon 
Systems  DSI-020/780  coprocessor  board  (4  Mbytes  of  RAM  and  equivalent  in 
speed  to  a  VAX-11/780)  in  about  8  hours. 
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